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1. INTRODUCTION 

Nowadays, the fifth-generation (5G) wireless networks is considered as a hot research area 
in academic and industrial based on the increasing data rate transmission. Relay communication, in which 
the relay forwards the signal received by a source to a destination, has a massive consideration in research, 
due to its ability to expand the coverage, increase the capacity, and reduce the power consumption [1-4]. 
Moreover, relaying the communication network is an effective way to combat the performance degradation 
caused by fading, shadowing, path loss, and is an efficient way to improve spectrum efficiency and extend 
coverage [1-6]. On another way, two-way relaying, where two users exchange information with each other via 
a single or multiple relays, provides improved spectral efficiency compared to conventional one-way relaying 
by using either superposition coding or physical layer network coding at relays. Furthermore, to satisfy the 5G 
requirements, relaying schemes with high spectrum efficiency, such as two-way, and full-duplex, etc., have 
been recently attracted considerable attention [7-10]. In [7], the authors investigate one-way full-duplex (FD) 
relaying and two-way half-duplex (HD) relaying to minimize the spectral efficiency loss associated with 
one-way HD relaying. In [8] and [9], the authors proposed and investigated one-way FD relaying with multiple 
antennas to provide a solution to overcome the spectral efficiency loss. In [10], the authors investigate 
one-way FD relaying with opportunistic relay selection to enhance the performance. 
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In this research, the two-way decode-and-forward (DF) HD relaying network with the direct link 
between two sources S; and S2 is proposed and investigated. The system model, energy harvesting (EH) 
and information transmission (IT) are presented in the first stage. After that, we derived the closed-form 
of the system outage probability (OP) in the next stage. Finally, the correctness of the analytical expressions 
is verified by Monte Carlo simulation in connection with the primary system parameters. The research results 
show that the analytical and simulation are the same in connection with the primary system parameters. 
In this research, we focus on some problems: 

a. The two-way DF HD relaying network with the direct link between two sources S; and S2 is proposed 
and investigated. 

b. The closed-form expression of the system OP is analyzed and derived. 

c. All the results are convinced by Monte Carlo simulation in connection with all primary 
system parameters. 

The structure of the rest of the paper can be drawn as follows. We provide the system model, 
the energy and information transfer phases in section 2. The closed-form expression of system OP analysis 
is derived in section 3. We introduce the results and have some discussions in section 4. In the last section, 
some conclusions are proposed. 


2. RESEARCH METHOD 

In this section, the two-way DF HD relaying network with the direct link between two sources S; and 
S2 is drawn in Figure 1. The energy harvesting (EH) and information transformation (IT) for this proposed 
model system are proposed in Figure 2. In this protocol, the transmission interval time is T, 
which consists of three-time slots T/3. In the first time slot T/3, the R harvests energy pP; from the source node 
Si, and the source uses the energy (1-p)P; for information transmission to R and Sz (here 0O<p<l 
is the power splitting factor). In the second interval time T/3, the R harvests energy pP2 from the source node 
S2, and the source S2 uses the energy (1-p)P2 for information transmission to R and S;. Finally, the remaining 
time slot T/3 is used for information transferring from the R to the nodes S; and S2 [11-16]. 





— EH Phase -—= IT Phase 


Figure 2. The energy harvesting (EH) 
and information transmission (IT) processes 


Figure 1. System model 


2.1. Energy harvesting phase 

Let S; transmits the symbol x; in the first interval time. We define the received signal at the R node 
R and the S node S: as: 

Vir = Nh 2X +n, 


I Aid I 
Vig =x, +n, 


d) 


Where Ella) }=R.E {>} is expectation operator and P; represents the average transmit power 


at the S;. Further, n! and n} denote the zero-mean additive white Gaussian noise (AWGN) with variance No. 
Then, the harvested energy at the R node can be formulated as the following: 
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In the second interval time, the source node Sz transmits the signal x2 to the nodes R and Sj. 
Therefore, the received signals at the R and S; can be expressed, respectively, as: 
IT IT 
Yor = Mo pXy tN, » 


P= I 
Yo = h, ,X; +n, 


(3) 


Where E{ln[} =F, P, is the average transmit power at the So. And we assume that n” and n/' 


are the zero-mean additive white Gaussian noise (AWGN) with variance No. Here, the total harvested energy 
at the R node can be defined as: 


2 2 
np(F|hra| + P,|hy -| )r 
Fe 


4 
h 3 (4) 


In this model, the average transmits power from source S; and S2 is assumed to equal. So, the (4) can be 
rewritten as the following: 


2 2 2 2 
_ne(P lA. + P, |h,,,| )r _ mpP (i.e +|h,,¢| )r 








5 
h 3 3 (5) 
Where we denote P, = P, = P.. Therefore, we have the average transmit power at the R as: 
E 2 2 
P, =—“*~ = P( +h 6 
8 (F/3) 1p. \*,.e| | sa i) 


2.2. Information transmission phase 
In the first interval time, after doing EH, S; broadcasts the information to the R node and S2 with 


remaining power (1—)P.. From that, we have the received signal at the R node and S> as: 


Vir = vl-p RX +n, 


; : (7) 
Viz = MX +N; 
Similar to the first case, the received signal at node and S; node can be given as the following: 
Yar = y1— phy 2X, Aa (8) 


Po Hl 
Yn, = My x, +N, 


Finally, in the third interval time, the received signal at the sources S; and S2 can be expressed, 
respectively, as: 


mn I 
y= he Xp +n, 


TI III (9) 
Vy = Ng Xp +14 
Where Ets} = P, . The received signal at the source node S; or S2 can be formulated as: 
yer = VJ Pehp ils @x)+tn ypr= VPrhp iC @x)+tn ypr= Jf Pehp i(x1 ® x2) + ny, (10) 
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Where i¢€ (1,2). Once again, after applying the canceling the self-interference method and combine with (6) 
the SNR at S; can be defined as: 


2 
h 
we Pal apy +{io,eP Joe =npvo\lor +4] 








(1) 
Where 9g, =|A,,| 
3. SYSTEM PERFORMANCE ANALYSIS 
Hence, the outage probability (OP) of DF mode at the S; can be formulated as: 
OP, = Pr{7n, <7, } = Pr{max {nove Lo, +9,].ve,} <7,3 (12) 
= Prinove, Le, +9,]<y,3xPriye, <7, } 


Where we denote: 


7 7 
th th 
Pr{ve, cag) Mth se oe 


3 


Which A, is mean of the random variable (RV) @;. Here we have: 


P, = Pr{npyeiesrth} = Pefon< Ze} 
nev 





J (14) 
= J Fy | d 
i on( ttt fo (pdp 


Where 9=9,+9, 


For calculating the probability in (13), the probability density function (PDF) 
and cumulative density function (CDF) of g is calculated as in Lemma | 
Lemma 1. 


The CDF of (9, + 2) can be formulated as: 


F(x) =Prlp< x] = Pr[g, < x-@Q,| 


x-% 


| £,@de, = {F, (x-Q,) fy, Qudy, 


0 0 


le X—-P, D4 
=— || 1-exp| -——— | Jexp| -— |d@, (5) 
ill [ A; ] [ i 
exo{-2] 
=1-exp — ee i) ee a) 
A, A We ae ; 


Where 4, is the mean of RV ¢,. 
Case 1: 2, =4, =A 
From (28), we have: 


xexp[-2] 
F(x) =1-exp{-%]}-. 42 


A 


= [f (Q,)dQ, 


(16) 
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And the PDF of @ can be computed as: 


x 
ress( =) 
Fe) A) (17) 


fo) ax 92 


Case 2: 2, # A, 


Similar, we can obtain CDF and PDF of g as followings respectively: 


a x |” “ * Vex x 
i z] Zale | z] : | «) oe 


worealalal| as) 


We can calculate the probability P3 into 2 case 
Case 1: 2, =4, =A 
Substituting (14) into (16), we can obtain: 








dy 
a? (20) 


1 co 
=1-- J ores Yth es{-2 0 
A“ 0 npyer, A 


Apply (3.471,9) in [17], (20) can be rewritten as: 


g 
orien 
Py =1- J exp 2} 4 


0 NPY PAL 











27th Vth 
=l- Ky] 2 21 
2 npyar, 2 npyar, oa 


Where K,(¢) is the modified Bessel function of the second kind and v"® order. 


Also, finally, substituting (13), (21) into (12), the system OP at the source S; in the first interval time can be 
obtained as: 


1 Vth 27th | Vth 
OP, 1 ea } 1- Kal 2 (22) 
BF BY npyar npyrd 


Case 2: 1, # A, 
Substituting (19) into (14), the probability P3 can be expressed as: 


oc 
Vth I @ 9 
Py i=1 [ex] } oa ca Jee 
3 0 newer | A4-A2 A4 Ay 
oc 
=l- : [ea “th } exo{ ©} eno{ 2] dp 
24-42 0 npyeA, 24 42 
oc 
1 : J ec “th pes ? Jew 
24-2 ¢ npyeA, 24 


1° Vth ? 
+ J e] 4 xexp dp 
24-/2 0 NPYVQA, Zo) 


Apply eq (3.324,1) in [17], (23) can be rewritten as: 




















(23) 
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npwAa, 
Pat fli = (24) 
—A2 \npy: 
ee Tojo [ig xK,| 2,| 2 
npwAag 


And then substituting (13), (24) into (12), the OP at the source S; in the second case can 
be obtained as: 


op2 | cxf - an )}- Z “th 
BY A4-AQ nev 





JFeve(2 Yth ) 


npyaAg 


= faxe(2 Vth ) 








(25) 





npwAa 


4. NUMERICAL RESULTS AND DISCUSSION 

In this section, we investigate the system OP using Monte Carlo simulation in connection with 
the primary system parameters [18-25]. In Figure 3, the effect of the energy coefficient 7 on the system OP 
is plotted with the primary system parameters as p=0.8, R=0.5 bps/Hz, y=0.5 dB. From the results, we can see 
that the OP of S; significantly decreases with the rising of n. Here, we consider both cases with 
and without a direct link between the sources S; and S2. From Figure 3, we can see that the system OP 
with the direct link is better than the system OP without the direct link between the sources. Furthermore, 
the simulation results agree well with the analytical results. In the same way, the system OP versus y 
is drawn in Figure 4. In this simulation, we set n=0.5, R=0.5bps/Hz, p=0.5, and consider both cases with 
and without a direct link between the sources. As shown in Figure 4, we can state that the system OP has 
a considerable decrease when yw varies from-10 to 10. In Figure 4, the simulation and analytical results 
are the same with all values ofa and b2, and again, the OP with the direct link is better than the cases without. 
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Figure 3. OP versus Figure 4. OP versus y 


Furthermore, we investigate the effect of R and p on the OP of the model system, as plotted in Figures 
5 and 6. In these Figures, we set the primary parameters, as shown in these Figures. From Figures 5 and 6, 
it can be observed that OP has a slight decrease when p varies from 0 to 1, and the system OP crucially increases 
with the rising of R from 0.5 to 5. In all two Figures, the analytical and simulation results agree well with each 
other. On another hand, the comparison of the system OP of both cases with and without a direct link between 
these sources is demonstrated in all Figures. From the results, we can state that the system OP with the direct 
link is better than another case. 
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Figure 5. OP versus R Figure 6. OP versus p 
CONCLUSION 


In this research, the two-way DF HD relaying network with the direct link between two sources S; 
S2 is proposed and investigated. The system model, EH and IT are presented in the first stage. 


After that, we derived the closed-form expression of the system OP in the next stage. Finally, the correctness 
of the analytical expressions is verified by Monte Carlo simulation in connection with the primary system 
parameters. The research results show that the analytical and simulation are the same as the primary system 
parameters. The results can provide a critical recommendation for the relaying communication network. 
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